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Abstract 

We propose a model of semi-vicious walkers, which interpolates be- 
tween the totally asymmetric simple exclusion process and the vicious 
walkers model, having the two as limiting cases. For this model we calcu- 
late the asymptotics of the survival probability for m particles and obtain 
a scaling function, which describes the transition from one limiting case 
to another. Then, we use a fluctuation-dissipation relation allowing us 
to reinterpret the result as the particle current generating function in the 
totally asymmetric simple exclusion process. Thus we obtain the particle 
current distribution asymptotically in the large time limit as the number 
of particles is fixed. The results apply to the large deviation scale as well 
as to the diffusive scale. In the latter we obtain a new universal distribu- 
tion, which has a skew non-Gaussian form. For m particles its asymptotic 

behavior is shown to be e 2^ as y ^ —00 and e 2™ y 2 as 
y — > (X. 

1 Introduction 

Exact solutions of l-dimensional (ID) many particle stochastic models [H] have 
given much insight into the physics of non-equiUbrium systems in one dimen- 
sion [35]. They serve as a testing ground for the macroscopic theories, being 
able to verify their predictions Examples are the description of different 
kinds of non-equihbrium phase transitions [TD], calculation of the large devia- 
tion functions for the density profile and total particle current [5], verification of 
the fluctuation dissipation relations [22 and testing of the range of their validity 

The range of models is very broad. In the context of the present article we 
mention two of them. The first one is the lock step model of vicious walkers 
(VW) that has been introduced in the physical literature by M. Fisher [2] 
to describe the wetting and melting phenomena. This is a random process 
defined as many non-interacting particles performing random walks on a ID 
lattice, whose space-time trajectories are forbidden to cross each other. The 
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term non-interacting means that the probabihty of a particular reahzation of 
the process, which meets the latter constraint, is given by the product of the 
probabilities of the random walks performed by each individual walker. The 
other realizations, where crossings occur, are assigned zero statistical weight. 
Such an elimination of a fraction of possible outcomes at every time step violates 
the probability conservation. A measure of the probability dissipation is the sum 
of the probabilities of all possible particle configurations at a given time, referred 
to as the survival probability. Its leading asymptotics for m particles has been 
shown by M. Fisher to decay with time t as a power law: t 

Another model, the totally asymmetric simple exclusion process (TASEP) 
PU] . has been widely discussed in connection to the Kardar-Parisi-Zhang uni- 
versality class [21] ■ In contrast to the VW model, this is a model of interacting 
random walks. The interaction prevents particles from jumping to occupied 
sites. Therefore, similarly to vicious walkers, the statistical ensemble includes 
only those events in which the space-time trajectories of particles do not cross. 
The difference is that there is an interaction that changes the statistical weights 
of particle trajectories when they pass via neighboring sites so that the total 
probability is conserved. In this case the quantity corresponding to the survival 
probability is just a probability normalization constant. 

Thus, the probability lost after imposing the global non-crossing constraint 
on the dynamics of non-interacting particles in VW is regained in the TASEP 
by adding the lacking probability locally at certain steps. In the present paper 
we consider the two models as limiting cases of a more general interaction, 
where the added probability is a varying parameter of the model that controls 
the probability dissipation, such that the probability conservation is restored 
when it is tuned to the TASEP value. In this connection a natural question 
arises: what happens with Fisher's asymptotics for the survival probability 
under such a generalization and, in particular, how does it cross over to the 
TASEP normalization constant. This is the first question we address in this 
paper. 

Specifically, we propose a semi-vicious walkers (SVW) model, which interpo- 
lates between VW and TASEP. It is a model of interacting particles with partial 
repulsion or attraction, where trajectory crossings are forbidden. The term par- 
tial repulsion (attraction) means that the probability for the particle to jump 
to an occupied site is not equal to zero like in TASEP but can be less (greater) 
than that of a free particle. At the same time, the non-crossing constraint leads 
to lack of probability conservation in the same way as in the VW model. The 
strength of the interaction, which also characterizes the probability dissipation, 
is a parameter of the model, which has the TASEP and VW as limiting cases 
at the endpoints of its range. 

In this article we obtain the large-time asymptotics of the survival proba- 
bility. Its limiting case corresponding to VW is given by the above mentioned 
result of Fisher, which yields the leading power law asymptotics. Later it was 
reproduced with more mathematical rigor together with the constant prefactor 
that was obtained for the particular initial configurations, where the particles 
are separated by equal spaces [25] . In the case of a general initial configuration 
of walkers this prefactor depends on the initial positions. This case has been 
studied in [55] . 

Our results can be roughly divided into two parts. For generic values of 
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the interaction strength, away from the point corresponding to the TASEP, the 
probabihty dissipation is finite. It is intuitively clear that the main asymptotics 
must be similar to the VW one. Indeed, we obtain the Fisher's power law with 
a constant prefactor that depends on the initial positions of the particles and 
on the interaction strength. It is shown to diverge in the TASEP limit. The 
second and probably the most interesting case is the transition region, which 
interpolates between the VW and TASEP behavior. To probe into this region, 
we consider a scaling limit of the survival probability, where the large time limit 
and the TASEP limit of the interaction strength are combined. In this way we 
obtain a scaling function of a single parameter that controls the transition from 
VW to TASEP. 

The second problem we address is the distribution of the integral particle 
current in TASEP. A first example of such an exact distribution has been ob- 
tained by Derrida and Lebowitz [B] , who found the large deviation function for 
the particle current in the TASEP confined to a ring. A specific property of 
the finite system is that there is a finite relaxation time, after which the system 
settles into a non-equilibrium stationary state, independent of initial conditions 
[IH] . Then, the tool used to study integral current fluctuations is, roughly speak- 
ing, an analysis of the relaxation of the system subject to a perturbation into 
the stationary state. Technically it is an analysis of the largest eigenvalue of the 
perturbed Markov matrix governing the process. 

In genuinely infinite systems the situation is more peculiar. In this case there 
is no characteristic relaxation time scale. When starting away from the station- 
ary state, the latter is never approached. In this case one needs to consider 
actual time evolution of quantities of interest, which depend on initial condi- 
tions. A major breakthrough in this direction has been achieved by Johansson, 
[2T] . He considered the TASEP evolution of an infinite cluster of particles, 
which initially occupies all sites of the lattice to the left of a fixed site, and 
calculated the distribution of the number of steps made by an arbitrary particle 
in this cluster. Johansson's solution has initiated a burst of activity in the field, 
which exploited deep connections of the TASEP to the theory of random matrix 
ensembles and the determinantal point processes. Results have been obtained 
for different initial conditions and extended to many particle joint distributions 
[31], [37], [53], [13], in,- Remarkably, in the scaling limit these results provide pa- 
rameter free universal distributions [35j of the fluctuations measured in the KPZ 
characteristic scale, which is of order of t^^^ as time t grows to infinity [26 . This 
is in contrast to the diffusive scale t^/'^, which, according to the Central Limit 
Theorem (CLT), characterizes the fluctuations of the distance travelled by a free 
particle [12j . The large deviation limit of the single particle current distribution 
has been studied in [19] in connection with the fluctuation dissipation relations. 

Despite the great success in finding the distributions of single particle cur- 
rents and their correlation functions, very few results on the integral particle 
current, i.e. on the distance travelled by all particles, are available for driven 
diffusive systems. In fact the only known exact result is the above mentioned 
large deviation function for the integral particle current for the TASEP in a 
ring [5] and its generalization for the partially asymmetric case [53], [H]. No 
results beyond the large deviation scale, neither a generalization for an infinite 
system has been proposed. On the other hand, extensive quantities like the in- 
tegral current, are important ingredients of the thermodynamics of the models. 
A knowledge of the character of their fluctuations could be of help for exten- 
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sion of the thcrmodynamical formalism to irreversible systems. The present 
paper makes a step in this direction. The problem we solve here is as follows. 
We study the large time asymptotics of the distribution of the total number of 
jumps made by a finite number of TASEP particles in an infinite lattice, given an 
arbitrary initial configuration. The idea that allows us to consider this problem 
in line with the previous one is the existence of a kind of fluctuation-dissipation 
relation that unifies the dissipation of probability in SVW and the statistics of 
fluctuations of the integrated particle current in TASEP. Specifically, an auxil- 
iary parameter, which violates the probability conservation, can be introduced 
into the evolution operator in TASEP to account for the total number of steps 
made by particles, see e.g. [8]. This parameter plays a role similar to the one 
played by the interaction in SVW, the two problems being equivalent after a 
certain change of variables. Using this fact, we interpret the result obtained for 
the survival probability in SVW as a generating function of the particle current 
in TASEP. The latter, in its turn, can be used to reconstruct the form of the 
current distribution. 

Like those for SVW, the results obtained for the TASEP particle current 
consist of two parts. The generic values of the interaction strength correspond to 
the distribution of the particle current at the large deviation scale, i.e. describes 
the deviations of order of time t. It turns out that it has a skew distribution 
with asymmetric negative and positive tails. These tails are connected by a 
middle part corresponding to the transition region. The latter yields the current 
distribution at the diffusive scale, i^/^, which is shown to have a skew non- 
Gaussian form, depending only on the total number of particles, and we suggest 
to be universal for particles performing a driven diffusion. 

One technical remark has to be made about the connection of our solution to 
the theory of random matrix ensembles. It is this connection which enabled the 
above mentioned progress in calculating the single particle current distributions 
and their many particle generalizations. In our solution this connection has 
also been exploited. Namely, the survival probability in the SVW model at 
generic values of the interaction strength and exactly at the TASEP point can 
be calculated in terms the Mehta integrals Im,k with k — 1/2 and fc = 1, 
which appear as normalization factors in the orthogonal and unitary Gaussian 
ensembles of random matrices respectively [55]. Note, however, that the scaling 
function obtained in the transition region for the system of m particles can be 
reduced to neither of these integrals except of at three limiting points, where 
it becomes Im. 1/2 Jm,i and Im-1,1 respectively. Thus, we obtain a generalized 
object, which interpolates between these three Mehta integrals, and, therefore, 
in a sense unifies three different matrix ensembles. To our knowledge no such 
generalization has appeared in the theory before. 

The article is organized as follows. In Section 2 we formulate the SVW 
model, state the results obtained and discuss their interpretation in terms of 
the probability distribution of the particle current in TASEP. Sections 3-5 are 
a technical part, where we prove the results outlined in Section 2. In Section 3 
we solve the master equation for the SVW model. In Section 4 we obtain the 
asymptotic formulas for the transition probabilities. In Section 5 we prove the 
limiting properties of the function characterizing the SVW to TASEP transition. 
Section 6 has a summary and conclusions. 
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2 Model and results 



2.1 Semi vicious walkers model 

Consider m particles on a ID infinite lattice. A configuration X of the system 
is specified by an m-tuple of strictly increasing integers 

X = {Xi < X2 < ■ ■ ■ < Xm} , (1) 

where Xi is the coordinate of i-th particle. The strictly increasing order reflects 

the exclusion condition, i.e. two particles cannot occupy the same site. We say 
that the relation X <Y holds for particle configurations if 

Xl < yi < X2 < ■ ■ ■ < Xm < Vm- (2) 

The SVW model is a random process, which is defined on a set of sequences of 
configurations X^, • • • , X*, such that 

X°<X^---<XK (3) 

We refer to such a sequence as a trajectory of the system traveled up to time t. 
Every such trajectory is realized with probability 

P(X°, ...,X') = T{X\X'~^) ■ ■ ■ T{X^, Xi)T(X\ X")Po(^")- (4) 

Po{X) is the initial probability of the configuration X and the transition prob- 
ability T{X,Y), from the configuration Y to X, is defined as follows 

rn— 1 

Y) = -d(xm - ym) ^ {xi - yi, Xi+i - yi) , (5) 
fe=i 

where 

i){k) = {l-p)6k,o+pSk,i, (6) 
e{k,l) = {l-p{l- KSi,i))6k,o+pSk,i, (7) 

and 

< p<l, (8) 
1-1/p < K < 1. (9) 

This means that at each discrete time step a particle can jump forward with 
probability p or stay put with probability 1 — p, provided that the next site is 
empty. If the next site is occupied, the probability for a particle to stay put 
is (1 — p{l — k)). The probability deficit p{l — k), corresponds to the process 
when the particle jumps to the adjoining occupied site, which is forbidden. This 
excluded process results in probability dissipation in this model. The form of the 
transition probabilities corresponds to the backward sequential update, i.e. the 
particles are updated starting from the rn-th particles one by one in backward 
direction. In particular limiting cases the model reduces to 

1. K = - VW, a particle jumps forward with probability p or stays with 
probability (1 — p), irrespective of whether the next site is occupied or 
not. But then those realizations of the process where two particles come 
to the same site must be removed from the statistical ensemble. 
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2. K = 1 - TASEP, a particle jumps forward with probability p or stays 
with probability (1 — p) provided the adjoining site is empty. When the 
next site is occupied the particle stays where it is with probability 1. 

The TASEP with the backward sequential update was studied in [2] and [37] , 
where it was referred to as a fragmentation model. In the case k = (1 — 
the probability for a particle to stay where it is when the next site occupied, is 
zero. Therefore, the trajectories of particles passing via neighboring sites have 
zero weight, i.e. they are removed from the ensemble as well as those which 
meet at the same site. Therefore this situation resembles the vicious walks of 
dimers. The range of k given in @ is due to the requirement for (1 — p(l — k)) to 
be a probability. Positive values of k correspond to repulsive interaction, while 
negative values correspond to an attractive interaction. The domain k > 1 is 
also of interest in connection with the current fluctuation in TASEP, though it 
does not have a probabilistic meaning in the context of SVW. 

2.2 The results about the SVW model 

Let us consider the quantity 

Vt{x'^)= P{X^...,X'), (10) 

x"<x^--<x* 

where the sum is over all the trajectories of the system starting at the configu- 
ration X^, i.e. Pq (X) = Sx,x°- This quantity is the partition function of the 
statistical ensemble of the trajectories with the statistical weights defined above. 
On the other hand, if we add the lacking processes allowing the particles to jump 
to an occupied site, the value of Vt {X'^) will have the meaning of probability 
for all the particles not to meet up to time t. In Fisher's original formulation 
of such a process, two particles annihilate when getting to the same site. Then, 
Vt (X^) is the probability for m particles to survive up to time t. Therefore, 
we refer to this quantity as a survival probability. Below we formulate three 
theorems, which specify the asymptotic behaviour of the survival probability in 
the limit of large time for different parts of the range of the parameter k. The 
proof of these theorems is the content of Sections 4,5. 

Remark 1 Two of the theorems below are stated and proved for complex val- 
ued parameter k and the third one for real k > 1. Obviously, the quantity 
obtained there has a meaning of the survival probability only for real k varying 
in the range Consideration of other values of k is justified by its later inter- 
pretation in terms of the generating function of the moments of the total particle 
current in the TASEP. In the latter case the complex values of k turn out to be 
meaningful and useful to reconstruct the total particle current distribution in 
the TASEP. 

2.2.1 The survival probability for SVW 

Generic case, |k| < 1 In this case the asymptotic behavior of the survival 
probability Vt {X°) as t — > oo is given by the following theorem. 
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Theorem 1 Let k G C 6e a fixed complex number from the domain \k\ < 1, 
and let — < oo for any i,j — 1, . . . ,m. Then, as t ^ oo the survival 
probability Vt {X^) for m particles is 

Vt{X^)^A{n-X'^)[tp{l-p)]-''^ [l + o({\ogtft-''^)] , (11) 
where the prefactor is given by 

A{.;X') = ^"nr^Ji/,^ + ^^ det [g.,{xl-x^;n)\ ^ (12) 
and where the function gij{x; k) is defined by 

5m(-; - f^^ ^ J, ■ (13) 

Thus, in the range k < 1, up to the factor A{X'^\k), which captures the 
dependence on the initial configuration X^ , the survival probability reproduces 
Fisher's power law. All the dependence on X'^ is in fact hidden in the deter- 
minantal part of A{X^]k). In some particular cases the determinant can be 
simplified to a more transparent expression. For example, for equidistant initial 
conditions, 

x'L-xl^aim-i), (14) 
where a is a positive integer, it can be calculated explicitly: 

det - a;°;K)]^^^ ^^^ = (a + K-oK)"*" (15) 

In the limit k — > the determinant reduces to 

detb,,(x^-x«;0)]= n (16) 

l<i<j<m 

Then, up to rescaling of space and time, one recovers the result [38j for VW: 



m m(m-2) 

7r^~2 J Y[ (2i~i)\ ' evenm 
A{Q;X')= n ■ ■ (17) 

^<'<i<'^ ^ 7r3-T2-^'^ n oddm 

In the limit k — s- 1, we have 

det[g,,,(x^-a;°;l)]^^^^^.^^ =1. (18) 

This limit corresponds to the TASEP. Hence the asymptotics must change, as 
the probability conservation is restored. The signature of this fact is the diver- 
gence of the term A{k; X°) that takes place in this limit. Specifically, (|12p and 
(fT5|) suggests that as n approaches one, A (k; X^) diverges as (1 — k)"™^™ 
Comparing the exponent of this expression with the one of the time decay 
^-m(m-i)/4^ we can guess that the transition takes place on the scale {1 — n) ^ 
1/Vt. This hypothesis is justified below. 
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Generic case, k > 1 In this case no values of k fall into the range 
Therefore, according to the Remark [1] the result formally obtained for Vt (X^) 
does not have a probabilistic meaning in terms of SVW. However, it is still 
meaningful for the description of current fluctuation in the TASEP. 

Theorem 2 Let kg (1,oo) be a fixed real number, and let ~ < oo for 
any i,j — 1,. . . ,m. Then, as t oo, the survival probability Vt for m 

particles is 



Here we use the common notations 



(19) 



l + O 



((logt)^t-i 



/2 



1 - 

i-q 



for the g-number and 



(20) 



(21) 



for the q- factorial. Note that the q- numbers turn to usual numbers [m]^ 
in the limit q — > 1. 

Transition regime k ^ 1 Consider the limit 

t — !■ oo, K — > 1, (1 — k) \/i = const. 
We introduce the parameter 



a ~ lim 

t — >oo 



^tp{l-p) 



(22) 



(23) 



. . , TO. 



Theorem 3 Let the condition \xl — Xj\ < oo hold for all i,j = 
Then,in the limit for the parameter a G C defined in 123\) taking any 

fixed complex value, Vt (^'^) converges to 



fm {a) 



1 + 



Vt 



(24) 



where the function (a) has the form of a multiple integral: 
fm (a) 



(-1)" 



(27r)~2!---(m-2)! 

oo 

dui J du2 

i=2 



dUrn / dV; 



dUr, 



(25) 



whe 



A{xi 



l<'i<j<m 



(26) 



is the Vandermonde determinant. 
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Figure 1: The log-log plot of the function f2{a) in the range a > 0. 



The argument of Vt in can be omitted as the dependence on the initial 
configuration is lost in the limit under consideration. The limiting behaviors of 
fm{oi) match the TASEP and VW asymptotics. Indeed, we prove in Section 5 
that 

lim a''^f^{a) = ——^\{T {1/2 + 1), (27) 

1=1 

U (0) = 1 (28) 
lim e-"^/™(a) = - -y (29) 

5Ra-+-oo [m — 1)! 

In (|27p and (|29p the imaginary part of a is implied to take an arbitrary fixed 
value. The proof of these three limits given in Section [5] is done by reducing 
fm (a) to different cases of Mehta integrals, Im.i/2, Im.,i and Im-i,i respectively, 
which are well known in the theory of Gaussian random matrix ensembles. |33| . 

The above results can be illustrated by the example of the two particle case, 
TO = 2, when the integral ([^5]) for /2 (a) simplifies significantly. 

/2(a) =e"'Erfc(a). (30) 

Here Erfc (a) is the complementary Error function 

00 

Erfc (a) = ^ / dxe"^'. (31) 

V'^ J 

a 

In Fig. [T] we show how /2 (a) interpolates between the SVW and TASEP limiting 
cases, which are (l?fl) and (1^ respectively. 



2.3 Current fluctuations in TASEP 

Consider the process with the transition weights T{X, Y) defined similarly to 
^ but where the functions d{k) and 6 (fc, I) are replaced by 

d{k) = (1-^)4,0 + 6-^^4.1, (32) 
eik,l) = (l-p(l-<5u))4,o + e^p4,i. (33) 
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Here < p < 1 and 7 is a complex- valued parameter. It is not difficult to see 
that these transition weights correspond to the TASEP except that for each 
particle, the probability p to jump is multiplied by an additional factor e''', i.e. 

e''''^PTASEp{X", . . . , X*) - f{X\ • • • f{X\X')Po{Xo), (34) 

where Ptasep{X° , ■ ■ ■ :X'^) is the probability for a sequence of particle config- 
urations . . . , X*, to occur in the TASEP for t successive steps and Yt is the 
total number of jumps made by all particles in this sequence of configurations. 
Thus, one can calculate the moment generating function for the cumulative 
particle current as follows, 

{^'"^')tasep= E e''^^PTASEPiX',...,X'). (35) 
x"<x^---<x* 

On the other hand, we can see that, if we define 

K = e^^, (36) 

P = n (37) 
(1 -p)e ^ + p 

then the following relation exists between the transition weights T(A, X') de- 
fined in (ISll),(ISni) and those of SVW, dH),©, 

(1 - p)-™ f{X, X') = (1 - p)-"' T(X, X'). (38) 

As a result we have 

{^'''')tasep = (1 - 1))"*^* m , (39) 

where Vt (A'°) is the survival probability calculated for the SVW model, and 
the parameters k, and p of SVW are related to the parameters p, 7 of TASEP by 
((S7)) . The function ((5^ encodes all information about the distribution of 
the integrated particle current. Thus, we can apply the Theorems to obtain 
the asymptotic form of this distribution. 

Large deviation function. It follows from the Theorems [T] and [H and the 
formula (f39|) that for fixed 7 e M the asymptotic form of the generating function 
of the particle current Yt is as follows 



TASEP 



[te^p(l-p)] i 



(1-^(1-6"^))*^^^ '—^ ; -T 7<0 

(40) 

From here we conclude that a scaled cumulant generating function of the random 
variable Ytjt exists 

A(7)= limllog(e^'-) ^ r mlog(l-fp(e--l)) 7>0 

^ - t-l™ i \^ I TASEP \ log (1 - p (1 - e"^)) 7 < • 
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It is convex and difFerentiable everywhere. Therefore, we refer to the Gartner- 
EUis theorem [TS],!!] to show that the random variable vt satisfies the large 
deviation principle with a rate function 

/ (v) = lim - logP {Yt/t = v) = sup (71; — A (7)) . 

The solution of the maximization problem yields 

J / \ f mB{v/m) V > mp 

\ B (v/m) V < mp ' 

where B {v) is the usual rate function of the Bernoulli process 

1 — p p 
B (v) =^ (1 - v) log h log - . 

1 — V V 

Central limit theorem scaling. The transition regime corresponds to the 
following scaling limit 

t 00, 7 — > 0, jVt ~ const. (41) 

To translate the results obtained for this case we consider the random variable 

y=lini^=^, (42) 
y'tpil-p) 

where the convergence is in distribution. Taking the limit (|1T|) in (j39p . using 
the Theorem [3] and noting that p —^ p as "/ —^ we obtain 

(e"^>TA5SP = e—^/^™ (a) , (43) 
where a is an arbitrary complex valued parameter related to 7 via - 



a = lim ^y/tp (1 - p). (44) 

t — ^00 

The random variable y is the rescaled deviation of the integrated current Yt from 
mpt, i.e. from the average value of Yt for m non-interacting particles jumping 
with probability p. Note that y is the variable that, in the case of free non- 
interacting particles, satisfies the conditions for the applicability of the Central 
Limit Theorem (CLT). According to the CLT the probability density function 
(PDF) of y for m independent particles is the Gauss distribution. 



P'„ny) = exp (-2/V (2m)) /V2^ (45) 
Correspondingly, the generating function of its moments is 

(e"^>frcc = exp {may 2) , (46) 

which is the first factor in the moment generating function Therefore, the 
form of the second factor, fmia), shows how the distribution of y differs from 
the one for free particles. 
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The moment generating function contains all the information about the orig- 
inal distribution. In particular, the cumulants of y are given by the derivatives 
of its logarithm at a = 0, 



{y")c-£^^og{e-y)^^sEPL=o- (47) 



The value of the first derivative, i.e., 

(2/>c = (j^>=/™(0), (48) 

shows how the difference between the mean velocity Vm of the center of mass of 
the particles and that of free non-interacting particles, which is p, decays with 
time t, 

mt y/t m 

A nonzero value of (0) implies that this difference is of order of t^^/"^ . As 
the TASEP interaction slows down the particle motion, one expects it to be 
negative, i.e. 

fL (0) < 0. (50) 

The second cumulant 

(y'), = (y') - {yf = m + (o) - {U^)f (5i) 

is related to the diffusion constant of the center of mass. 



1 {Y?) {Y,y p{i-p) I 
^™ = — ^ — = — 2 — \y )c (52) 

The next cumulants, e.g. 

{y^)^ ^ (y3)„3(y^)(y)+2(y)^^ (log/™(a))"'|_o, (53) 
{y')^ ^ {y')~4{y'){y)-3{y')' (54) 
+12 (y)^ (2/2)- 6 (log /„(«))(*) 



a=0 



quantify the discrepancy of the distribution from a Gaussian form, being iden- 
tically zero for the latter. 

The asymptotical behavior of the generating function at large absolute values 
of Sfta can be readily obtained from the ones of /„(«), ([27|) . ([29|) . 



- /TASEP 



a 



The PDF of the random variable y can be obtained as an inverse Laplace trans- 
form of its moment generating function (|43p 



/3+ioo 

Prniy)^ J e™"'/2-"V™ (a) ^ (56) 

/3— ioo 
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As the function /,„ (a) is bounded and analytic in any vertical strip of finite 
width, the parameter [3 can be chosen arbitrarily. The asymptotic results (j55p 
for the generating function can be used in the integral ((5^ to evaluate the 



asymptotics for PDF Pm{y)- Choosing (3 — y / m ioi y 
y ~f —oo we obtain 



oo and (3 — y/r 



for 



. n"ir(i/2+i) 2" 



(57) 



/27r(m-l)! 



V 



Thus, the form of the distribution Pm{y) is far from being symmetric, having 
tails of two Gaussian-like functions with different dispersions, rn? and m, on 
the left and right respectively, the latter also multiplied by "Fisher's factor" 

y — m(m — 1)/2 

Let us compare these result with the data obtained from Monte Carlo simu- 
lations. We modelled the TASEP for to = 2, 3, 4, 5 particles, which have evolved 
for t ~ 10^ time-steps, the statistics having been collected from 10^ samples. We 
would like to compare the data obtained for the generating function (e"^) and 
the PDF Pm {y) with our predictions. An explicit evaluation of these functions 
requires detailed analysis of the function fm (a), which is given by the multiple 
integral For arbitrary to this needs a significant calculational effort, which 
is beyond the goals of the present article. Fortunately, for to = 2 the function 
/2 (a) is simple enough, being given by pop . and we can use it for plotting 
the generating function and the distribution. In Fig. [2] we show a plot of the 
logarithm of the m — 2 moment generating function, whose analytic expression 
is 

(e"'>TASi;p = e2"'(l-Erf(a)). (58) 

It has a skew convex form growing more rapidly to the left than to the right, with 
a minimum at a = 0.432752. One can see good agreement with the numerical 
data in the central part of the graph. There is some discrepancy at the tails, 
which can be attributed to the finite-time corrections, i.e. the lack of statistics 
of large events at the finite period of measurement, which becomes significant 
when the absolute value of a is large. 

The function (|58|) allows a calculation of any derivatives, and, hence, of any 
cumulants of the random variable y. In Table we show the first four cumu- 
lants for to = 2, the case of two particles. Their values are in good agreement 
with the results from the Monte Carlo simulations. In Fig. [3] we show the result 





Analytic 


Numerical Analytic 


Monte Carlo 


{y)c 


-2^-1/2 


-1.12838 


-1.12545 


{y'l 


4 - 47r-i 


2.72676 


2.72518 


(y')r 


4(7r - 4)7r-3/2 


-0.616636 


-0.617642 




32(7r- 3)7r-2 


0.459083 


0.498263 



Table 1: Cumulants of the random variable y. 



of numerical evaluation of the integral (|56p for to = 2. There is a very good 
agreement with the simulation results. At first glance the form of the distribu- 
tion shown on Fig.[3K appears Gaussian-like. A more accurate impression of the 
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Figure 2: Plot of the logarithm of the moment generating function for m = 2. 
Solid line is the plot of the formula ([58]) . Dotted line is the result of Monte 
Carlo simulations. 




Figure 3: Probability density function Pm{y) for m = 2 particles (a) and its 
logarithm (b). The solid line shows the theoretical predictions for the distribu- 
tions. 

form of the distribution is given by the logarithmic plot of Fig. ^jp which shows 
that the distribution is actually skew, decreasing more rapidly as y grows than 
as it decreases. 

Simulation results obtained for more than two particles can be tested against 
the asymptotical formulas (l57|) for the tails of the distribution Pm (y). In Fig. [4] 
we plot the distributions measured for to = 2, 3, 4, 5 particles, (Fig. 2^), and its 
logarithm, (Fig. [Hd), the latter being compared with the graphs of ([S7| . 

One can see that for all the four graphs the left tails are perfectly fitted 
already for rather small values of y. A good fit of the right tail takes place only 
for TO = 2. For to = 3 only a few data points approach the asymptotical line, 
i.e. in this case the right asymptotical regime is actually at the borderline of 
the statistics. In the cases m = 3, 4 the statistics available is clearly not good 
enough to reach the asymptotical regime. Significantly larger evolution time 
and statistics would have to be considered. 
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Figure 4: Probability density function Pm{y) for m = 2, 3, 4, 5 particles (a) and 
its logarithm (b). The solid lines show the theoretical asymptotics of the tails 
of the distributions. 

3 The master equation 

From now on we consider only the SVW model dependent on the parameters 
p and K. Our first step is to calculate the probability Pt {X, X^^ of transition 
from the configuration to X for arbitrary time t: 

Pt{X,X")= P{X\...,X'). (59) 

X0<X1---<X*eeX 

The method of finding the transition probability was first developed by Schiitz 
for the continuous time TASEP ^313], who used the Bethc Ansatz first applied 
to the ASEP by Gwa and Spohn 18J. Here we follow a similar procedure. The 
transition probability obeys the master equation 

Pt{X,X')^Y.'^{X,X')Pt.,{X',X% (60) 

X' 

the transition weights T (X, X') being defined as above, ©-(HI). The problem of 
finding the eigenvectors and eigenvalues of the matrix T (AT, X') can be solved 
by the Bethe Ansatz technique. As this technique is rather standard and has 
been reviewed in many monographs, we simply state the results here. For details 
of similar derivations, the reader can consult for example with the review [5]. 
As a result we obtain the solution of the left and right eigenvalue problems for 
the Markov matrix T {X, X'): 

K{Z)^>z{X) = ^T(X,A')*z(A'), (61) 

X' 

K{Z)^z{X) = Y.T{X',X)^z{X') (62) 

X' 

parametrized by an m-tuple of complex parameters Z = {zi, • • • , z„i} . The 
corresponding eigenvalue is expressed in terms of these parameters, 

m 

A(Z) = [](l~p + p/z,), (63) 
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and the eigenvectors are given by the foUowing determinants 

^z{X) = det(zr^(l-.^z,)^-^)i<,,<„,, 
'^z{X) = det(2p^(l-Kz,)^- 



l<'i,j <m 



(64) 
(65) 



It is not difhcuh to check that these two eigenfunctions can be used to construct 
the resolution of the identity operator 



27riz,- 



= S 



X,X' 



(66) 



where the integration over each z.^, « = 1, . . . , m, is along a contour of integration 
which has to satisfy the requirement that the pole of the wave function at 
z = 1/k has to lie in the exterior. Then the solution of the initial value problem 
for the master equation is given by 



Pt {X, X")^—J A* (Z) {X) ^z n ■ 

i—l 



27:1 Zi 



(67) 



Finally we end up with the following integral expression for the transition prob- 
ability 



A*(z)n 



X det z. 



1=1 

^0 ^0 



{l~nz,f-^ 

(1 - >^z,y-^ 



(68) 



n 



l<i,j<m i—\ 



dz, 
2'Kizi 



The integration can be easily performed by counting the residues. The result is 
a determinant of an m x m matrix of the form similar to the one obtained for 
the discrete time TASEP with backward update ( [2] , [37]). Note that in the 
case of vicious walkers, k = 0, the eigenfunctions are of free fermion type 



-^ziX) = det(zf^), . , , 
'^z{X) = detfzr"^) 



(69) 
(70) 



and the integration yields the famous Lindstrom-Gessel-Viennot theorem [50] . 

P, {X, X°) = det [Fo (x, - t)] ^^^^^^^ , (71) 



m- 

where 



Fo{x,t)=p^l-pf 



(72) 



These formulas serve as a starting point for the asymptotical analysis of the 
survival probability. 
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4 Asymptotic form of the survival probability 



To obtain the survival probability Vt {X^) for SVW we have to sum the tran- 
sition probability Pt (X, X") over the set of all final configurations X: 

(73) 

{X} 

We solve this problem in the limit t ^ oo. For pedagogical reasons we first out- 
line the derivation for the VW model, which simply reproduce Rubey's results, 
|38j . The procedure we use amounts to an asymptotical analysis of the expres- 
sion for Pt (X, X^^ by means of the saddle point approximation for the integral 
()67p . which reduces the sum over final configurations to known integrals. The 
main ingredients of the derivation for the VW model are then applied similarly 
to the SVW model but with some modifications. 



4.1 Vicious walkers 

In the case of VW (k = 0), the integral ([SS]) takes the form 



/ A* (Z) n 4'"^- <^{4--%<..<m fl (74) 
fc=i 1=1 '■ 

Here, the determinant under the integral can be expressed in terms of the Van- 
dermonde determinant 

A(Z) = det(zr"')i<»j<m= n (75) 

l<'i<j<m 

and the Schur function 

s^{z,, . . . , z™) = det (zf'+"-') /A {Z) (76) 
parametrized by the partition x = (xi ^ ' ' ' ^ Xrn > 0) defined by 

X = (4. - re? - m + 1, a;?„ - - m + 2, • • •), (77) 

as follows: ^ ^ 

det(zf )i<»j<™ = A (Z) {Z) . (78) 

Thus ((74|) can be rewritten in the following form 

PT{X^X^)=U{Z)s^{Z)X{e'^'^^^'^f^, (79) 

•' i=i ' 

where 

h,{z) = \og{l-p + p/z)+v,\ogz (80) 

and 



V, = (81) 
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Now we are ready to estimate the integral asymptotically as i — > oo. We assmne 
that the differences (a;° — a:°) are kept bounded for any i and j. The saddle 
point of the function under the integral is defined by the equation 



which yields 



Zi = 



(1 - Vi)P 



(1 -p)vi' 

In the vicinity of the saddle point hi(z) has an expansion 



h^{z*+0 = log 



1-P 
1 - Vi 



1-Vi 



(82) 



(83) 



(84) 



4 C-^Xj^e+oie). 

2 V P ) 1 - "" i 

The integration contours can be deformed to a circle centered at 0, crossing the 
real axis at z*. Writing points on the circle as Zi = 2:*e'"^', we have 



^{h,{zC)) - h{z*) + i log [(1 - v,{\ - cos((/.,))2 + V. 



2sin2( 



0] ■ (85) 



It follows that there is a single maximum at = 0. Moreover, since all deriva- 
tives are bounded provided Vi < 0, the saddle-point approximation holds uni- 
formly in Xi. (Note that ((74|) is zero if > 1, and if Dm = 1 then the probability 



p can be extracted as a factor, the remaining integral over zi, 



being 



of the same form.) It is easy to see that the contribution to the sum over X 
from points with Vi > p + e for any fixed e > is negligible in the limit t ^ oo. 
(Below, we shall see that the effective range of the summation is in fact even 
smaller.) The saddle-point approximation [11 now yields 



Pt(X,X°) = 



n 

1<? <j<m 



i-p 



— ti "e ^ (i"p>" 



n 



(1 " Vx)p (1 - Vm)p 

{\-p)vx""' (\-p)Vm 



(86) 



l + O 



The next step is to perform the summation (|73p over the range of the final 
configurations X S \x\ < x\ < ■ ■ ■ < Xm < oo}. For this we need to demon- 
strate that holds uniformly in X . To this end we first show that the main 
contribution to the sum comes from the domain 



pt — Vilogt < xi < ■ ■ ■ < Xrn < pt + Vilogt. 



(87) 



Indeed, hi (^j^z^^^ is a concave function of Vi in the domain Vi £ (0, 1) with a 
single maximum Vi = p. It follows then for — pij > vilogt 



thi 



< e 



th ( i-^'°g'/(i-P) ^ 
'V l+Vtiogt/p J 



e 2p(i-p) 



l + O 



logi 
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All the other factors in ([55]) are at most of polynomial order in t, while the total 
number of nonzero terms in the sum of interest ((7S)l is O (i™). Therefore, the 

contribution from the complement of (|87| being of order of O I t^e 2p(i-p) i fQj. 

some constant s is asymptotically negligible compared to the contribution from 

In the latter one can approximate the function hi ^ [i^p')^^ ^ by the second 
term of its Taylor expansion at vi — p, which yields 



= ^TF^x (1, • ■ • , 1) 

(2^)f {tp{l~p)) - 



(89) 



J|exp 

i=l 



(xi - x'^^-pty 



n ~ 

l<'i<j<m 



1 + 



2ip(l-rt 

We now have to evaluate following sum in the limit t — > oo, 



J2 Yle'^"^ n (90) 

-Vtlogt<a:i<---<2:,„< \/tlogt*=l l<j<j<m 

This can be done by means of the following lemma: 

Lemma 1 Let h : M™ R be a twice differ entiahle function of at most poly- 
nomial growth. Then as S 0, 



5^ Myi,---,2/«)n^""'' (91) 

yi<...<ym.\yi€SZ t^l 

dyi / dy2... dy^h{yi, y,„)e"5 i vl _,_ o{S). 

-OO Jyi Jym-l 

Proof We subdivide the domain — oo < xi < . . . < x^. < oo into hypercubes 
of the form 



Bsiyi, ■ ■ ■,yrn) = {(a;i, . . . ,a;m) : max|a;i - y,] < -}, 



(92) 



where j/i € (5Z and yi < ?/2 < ■ • ■ < Um- The remaining region is small and 
its contribution will be estimated shortly. We then write, for (xi, . . . ,Xm) G 
Bsiyi, ■ ■ • ,ym), 



hixu x™)e-^(-?+-+-") - h{y,, . . . , y„)e-ife?+-+«™) 

d 



(93) 



< sup 

(ui,...,ti„)eSa(j/i,...,j/m) 



max 

i=l 



du. 



Mui,...,u^)e-i('^'+-+<^ 



\x^ - yi\ 



Now, 



dU; 



■/i(wi,...,"m)e-5("'+-+"™) 



(94) 



(^-^h{ui, . . . ,Um) - Uih{ui, . . . ,Um)^ 6 
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which is easily seen to be bounded by C'e^^^i+ ' +^"i)/2 for some constant C > 0. 
It foUows from the convergence of the sum X^yssz '^^"^ uniformly in S that 
the difference between the integral over the region 



U Bsivi, 

Vi<...<ym\Vi&S1' 



,Vni) 



(95) 



and the sum is of order 5. There remains the integral over the complementary 
region, but this is obviously of order 5 as the integral converges and the region 
has width 5. 



Then, after going to rescaled variables yi — Xi/ \/tp{\ — p) and writing 5 
1 1 \Jt-p{\ — p) the sum (j90| reduces to the integral 



{tp{\-p)) - 



dXr, 



n 



l<'i<j<m 



(96) 

(the range of summation is extended to {~oo < xi < ■ ■ ■ < Xm < oo) by the 
same argument as above. Note that the absolute value signs \xj — Xi\, though 
redundant in this range, are nevertheless useful as they make the expression 
symmetric with respect to permutations of the variables xi, . . . , Xm- One, then, 
can use this fact to extend the integration to the whole M™, which yields an 
additional factor of m!, which has to be compensated in the end. As a result 
we obtain 



1 



,1/2 



where 



[p{l-p)t\ 



dyp ■ 



dy2 / rfj/iexp 



1) 



l + O 



Vi 



1 Y^'" 2 

2 z^,:=i 



(97) 



X n 

l<i<j<~nL 



,2k 



m/2 -fr r(/fc+l) 



(2.)™/^n 



1=1 



r(A: + l) 



is the Mehta integral [33], which first appeared in the context of Gaussian ran- 
dom matrix ensembles. Finally, one can use the following formula for the Schur 
function [31] 



s^(l,...,l) 



n 

l<i<_7<m 



resulting in the following expression for the survival probability: 



Vt (X°) 



1 



n 



r (1/2 + 1) 



b(i-p)t]™(™~')/^^"/' ^^-^ ^! 

Vi 



(99) 



(100) 



X n ( 

l<i<j<rn 



Xj X^ ) 



1 + 
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After reexpressing the gamma functions in terms of factorials we obtain the 
form given in (|17p . 



4.2 Semi-vicious walkers 
4.2.1 The case of generic 

To study the asymptotic behaviour of the survival probability for the case of 
general k, one can start with the following integral representation for the tran- 
sition probability 



i=l 
X 



X° -X° + l Xi-X° 
Z; 



^*(^) n n (101) 

l<ij<rn ^ l<i<j<m 

where the integration in each variable is along a small circle around zero, \zi\ < 

\ for any i,j = l,...,m. This representation can be reduced to the form 
by direct integration over each S,j (j = 1, . . . , m). This is done by summing the 
contributions to the integral coming from all the poles = Zi, i — 1. . . . , m. 

Though the most of analysis of the large t asymptotics of this expression 
is similar to the one for VW, one important difference exists. The expressions 
under the integrals over z^, i = 2, . . . ,to, have singularities at Zi — 1/k, the 
poles of the form (1 — KZi)^~^ , which can be located between the origin and the 
saddle point. In this case the contour being deformed to the steepest descent 
one, crosses this singularity and its contribution must then be extracted from the 
saddle point contribution. While for |k| < 1 this does not affect the asymptotics 
of the sum over xi, . . . , Xm, evaluated subsequently, for |k| > 1 its contribution 
turns out to be dominant. 

It is, however, technically difficult to calculate the residue at the multiple 
pole of the complicated expression. To avoid this calculation and to evaluate 
both cases in one go, we expand the term (1 — KZi)^ ' into a series in powers of 
(KZi) and then integrate it term by term in the saddle point approximation. As 
a result we obtain Pt (X, in the form of an (m — 1) - fold series 

fe=i -'Co 27ria 

l<i<j<m {n2,...,n„}eZ!J-i 



where 



11 /o , n ^ 11 '^''^-''^^ 
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X 

l<z,j<m 



^1 = (104) 



The next step is to use this approximation to perform the summation of p02p 
over the domain {x^ < xi < X2 < ■ ■ ■ < Xm < 00} . The effective range of this 
summation depends cruciaUy on the behaviour of the other sum in n2, . . . , rim- 
Namely, the effective summation range is different depending on whether the 
value of K is greater or less than one, when the term k"' is decreasing or in- 
creasing respectively. We consider these two cases separately. 

The case |k| < 1 Here k takes arbitrary complex values in the domain 
|k| < 1. As in the case of vicious walkers, we argue that the exponential part 

makes the whole expression negligible beyond the range 
p - log i < < p + log t. (105) 



exp 



th, 



In addition, for Ui > (logi)^ / |log |k|| we have 

k"' ^ =Oit-'°s*{\ogtf^'-^'>), (106) 

so that we can limit the summation over rii to rii < (logi)^ / |log|K||. There- 
fore, Ui is negligible compared to Xi in the domain (llOSp and can be ne- 
glected in the definition (|104p of Vi. In the range (|105p we can approximate 
-4 ({^i, Xi/t}^^ , {?ifc}™ 2 ' ) tiy the leading term of its Taylor expansion at Vi = p, 
i ~ 1, . . . ,m. 

-4({6,a;»A}"i >{"fc}™2) = 

' ^ u^- '^:^~') (107) 



tp{l-p)j fj^ \ n 



^=l V2^ip(l-P) i<,<,<„, 



X 

i=l 



\\ogtf 



One can see that the terms dependent on {xi} and {tii} decouple and the terms 
dependent on 712, ... , rim can be summed up. 
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The remaining sum over Xi for i — 1, . . . , to is transformed to an integral using 
Lemma [TJ 

^ ^ A{{i.,x,/t}Z,Ank}T=2) (109) 

x1<xi<X2<---<x,n {n2,---,nm}eZ™Q 



x-1 



{tpil-p)) (27r) 



11(^^-1) 



2^1 



i + o 



Combining dZS]), (fT02ll and (fT09| we obtain 



'ma/2 



^6 



n 



jr> 27rif, V 1 — K 



m y-x"-x" + l 



1 + 



n fe -^^ 

l<?<j<rn 



Vt 



(110) 



where Im,i/2 is the Mehta integral defined in Writing the above product 

of integrals in determinant form and using the definition of Im,i/2 we arrive at 
the final result 



Vt {x°) 



2™ _ 

b(i - p)t 



2/2 



n r (V2 + 1) det [(g,,, (4, - 4))" = J • (111) 



1=1 



Here the function gij{x) is defined as follows 



{n + n^- (1 + 



Co 



27ri 



(112) 



The case k > 1 Let k be a real number, k > 1. We return to the formulas 
(|102|103|) . The crucial distinction from the case |k| < 1 is that the presence of 
exponentially growing terms k"* affects the range of values of wi , . . . , Vm , which 
make the major contribution to the final sum of p02p . Indeed, we can write 

k"' = e*^''°e''K-^'+='". (113) 

Therefore, if we keep Vi fixed, the sum over Xi is rapidly converging. At the same 
time the maximum of the Wi-dependent exponential part of the r.h.s. of (|102p 
is shifted due to the appearance of the additional term tvi log k. In a sense, the 
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roles of the variables X2, . . . , Xm and n2, . . . , rim are interchanged compared to 
the case k < 1. 

Consequently, instead of summing over 712, ... , and then over . . . , Xm 
we go to the variables vi, . . . ,Vm, (|104p . and X2, ■ ■ ■ , Xm and evaluate the sum 
over the latter first. 

E E ^ ({6,^ Jill, w}™ 2) (114) 



E 

{ill, 



V2t+X 



Vmt + X" 



E ••• E 

2:2— a^l + l Xm—Xj-n-i-]-! 



A 



Collecting the factors of A {{£.i,Vi}^i , {vi — {xi — a;°„) /t}^^) dependent on 
X2, ■ ■ ■ , Xm we can evaluate the sum over these variables 



V2t + xl^ 

E 



Vmt + X" 



E n- 



-Xi+x" 



X2—X\+l Xjn—X„ 



1 + 1 i=2 



^(™-l)(.l-xi) ^ 



(k - 1) • • • (k*^ 



)"('-2)! 



tvi — Xi -\- -h i — 2 

tVi Xi "|~ X™ 



(115) 



Here we extended the upper limit of all the summations to infinity, which yields 
a correction of order of and we used Stirling's formula to approximate 

the binomial coefficient 



I + n 
n 



1 + 



(116) 



We also imply that the value of Vi in the effective summation range is finite and 
positive. Indeed, the range of summation over Vi is defined as above by the re- 
quirement that the exponential parts of A {{£,i,Vi}^i , {vi — {xi — a;°„) /Oi-2) 
are not too small. Specifically, the exponentiated expressions are 



exp < t 



for i = 2, . . . , m, and 



exp < t 



hi 



(1 - Vi)p 

{1-p) Vi 

(1 - Vi)p ' 
(1 - P) Vi 



+ Vi log K 



(to — 1) Vi log K 



(117) 



(118) 



the term t{l — to) vi log k — log ^k'™ i)(^m ^1)^ jjj ^jjg latter coming from the 

result of the summation over X2, ■ ■ ■ , Xm, (|115p . For a real positive k the major 
contribution to the sums over vi, . . . ,Vm comes from the neighborhood of the 
maxima of the exponentiated expressions 



\Vi - Ui 



< log* 



(119) 
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where the maxima Ui are located at 

ui = u(k1-™) (120) 



and 

for i = 2, . . . , TO where 



= u(k) (121) 

Then we can follow the above procedure to evaluate the sums over vi, . . . , Vm- 
Going from the sum to an integral over the variables 

y,^V-t^^=^. (123) 
l/Ui (1 - Ui) 

we arrive at the integral expression 

E E ••• E 

>X J 



X 



(u — Ul)™ (1 — u) 



= (1-p + kp)^"-'^^(1-p + k1-»^ 

(m-l)(m-2) rn(m-l) 



X ■ 



ym(m-l)/2y™-l \l-p 



-\-OQ +00 



X 

5-2 



2<i<j<m 



ni»-;^ (-«(t))- 



Note that we keep the leading terms of the Taylor expansion in {vi — Ui), i = 
1, . . . , TO, everywhere under the integral except the product in the last line, where 
we keep the terms of two subsequent orders. The reason for the latter is that the 
leading terms of the multipliers cancel due to the antisymmetry of the rest of 
the expression in W2, . . . ,fm. Therefore, what contributes is the antisymmetric 
part of this line, that is Y[2<i<j<m ivi ~ Vj) / ('^ ~ !)■' which contains only the 
terms of the same order. Inserting it, we again arrive at the Mehta integral 
Im-1,1 over (to — 1), variables 2/2, •■ • ,2/m, while the integral over yi decouple 
being just the Laplace integral. After substitution of the explicit form of u and 
ui the r.h.s. of p24p becomes 



{1-P+ Kpy"'-'>' {1-P+ K^-"'p) (k" - 1)™"' K 
m 

n [(*- 1)1(^-1-1)] 

1=2 
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m / ^\ 1— m 

,i(2^)-'^nUfc"T (125) 



1— m 

This formula together with ((731), ITM)) and ([M]) yields 



i=l 



n fe-^^)- 



(e, -«,-!)— (e,-^"-i) 

The integration over ^fc, fc = 1, . . . , m is performed over a circle encircling the 
singularities of the expression under the integral, i.e. |^fe| > k™"^. First we 
integrate over then over ■ • • , ^i- It turns out that, if we integrate 

in this order, the expression under the integral being evaluated will contain only 
one simple pole at each step. As a result we arrive at the simple one term 
expression 



TT ^ — TT (c,-c.)(i27) 

= (-i)^^«,i:rj-/(-s?-°J 

This formula together with (1126p results in . 

Remark 2 An extension of the approximation technique used to complex val- 
ues of K is problematic for |k| > 1. The reason is that in this case the critical 
points ui, . . . , Ura are away from the real axis. It follows then that their contri- 
bution can be exponentially smaller the other corrections appearing. 

4.2.2 The limiting case k — > 1 

Now we consider the limiting case 

t —> cxD, K — ^ 1, (1 — k) ^/i = const. (128) 

We start with the formula and expand the term (1 — KZi)~'+^ into its Taylor 
series, 

\ / l<ij<m 

The integral in the determinant can be evaluated in the saddle point approx- 
imation, the analysis being similar to the one above, ([5n |) - ([55|) . with the same 
function hi{z), ([80|l except that Vi now depends on nf. 



(130) 
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What is special about the hmit k ^ 1 is that the saddle point can coincide with 
a zero of the factor (1 — nzy within the effective range of the summation over 
Vi. Therefore, instead of expanding this term into a Taylor series, we leave it 
in the integral as is, while the rest can be expanded around the saddle point 
as usual. Then we use the following formula for Hcrmite polynomials (see [17j . 
formula 3.462.4), 

OO ^ 

J e--" {x - (3f dx = (^^^ i?n(i/3). (131) 



As a result we obtain 
dz 



JhA^)z-".-"(l - nzy-' = e^'^^^'Uzy-'^^-' (132) 

ZTTIZ 



r ^e-5*l""(^')l«'(l-Aj(z*+iO)-''~' 



thi 



P 



1-pJ ^a,?„-.0-3j72 



l + O 




Next, we argue that the dominant range of the summation over X and 
is the domain 

pt-Vi\ogt< Xi + rii <pt + Vilogt, (133) 
where Xi varies within the range 

Xi <xi < ■ ■ ■ <x,n <t, (134) 

and 

< < oo. (135) 

To this end, consider the integral p32p for some particular i and j. After 
expanding (1 — KZi)-' into a binomial sum, it becomes a finite sum of terms like 



^^~P> P [ n, + x,~x°^+k ) 



(136) 



where fc is a finite integer, < fc < j. Beyond the range (|133p this can be 

. _ (log i)^ 

estimated using Stirling's formula to be 0(t~ ' e 2p(i-p))^ xhe summation 
over Xi, which includes at most t nonzero terms, multiplies this estimate by a 
factor of t. Finally the summation over rii yields an additional factor (1 — 
with the result that the order of the contribution from outside the domain (|133|) 
is 

o(^il-Kr't'^^e'^^y (137) 

Below, the leading term of the sum of interest will be shown to decay at most 
as a power of t. Therefore, when k is such that (1 — k) = ©(i"") with any fixed 
s > 0, the term (|137p is asymptotically negligible. 
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One can approximate (|132p using the Taylor formula, which yields 



det 



ltp{l -~p) f Vi-p 



p{l-p) 




TT^ [2tp {p - 1)] 



l + O 



Vi 



(138) 



The determinant can be simplified by adding to every line a multiple of the lines 
below it, such that all the terms of the Hermite polynomials except the highest 
cancel: 

det (aO]i;;-=i = (-2)'^'^ A(ai, . . . ,a^) (139) 

Thus, the survival probability takes the following form 



(-1)- 



E E — _ 

{-.irii 2 (2^) ' [p (p - 1)]" 

™ _ (..+..-.;k-.0^ f ^ + n,-2 

i=2 



(140) 



X A{xi,X2+n2,...,XrrL+nm) 



l + O 



where the summation is over the domains of {a;;}™ and {n^jflj defined by the 
inequalities (|133p - (|135p . Due to the presence of the Gaussian factor exp{—t{vi — 
p)^/(2p(l — p)t)) the sum over rii converges uniformly in Xi. Therefore we can 
interchange the order of summations over Xi and over rii. This allows us to 
apply Lemma [T] first to the variables Xi and then also to the variables n^. As 
the characteristic scale of Ui is of order ^/i we can use the approximation (1116p 
for the binomial coefficient, where the correction term yields an error of the 
order O (i~^/^) in the final result. To write down the final limiting formula as 
t ^ oo for P (X'^) we introduce the rescaled variables 



Ui 



(Xj ~ pt) 

^/tp{l~p)' 
rij 

Vtp{l-p)' 



(141) 
(142) 



and the transition parameter a, (|23p . which is constant in the limit under con- 
sideration. The formula ()140p then takes the form , where (a) is given 
by 125]). 



5 Asymptotic behaviour of fm (a) 

In this section we evaluate the limiting behaviour of fm{a) for a — > cxd and 
a — > — oo and its value at a = . In the latter case it is just the probability 
normalization of the TASEP, so fm (0) must be equal 1. The limit a — > — oo 
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has no probabilistic meaning, but it can be considered a particular limit of the 
generating function of the rescaled particle current in the TASEP: see Section 
2. Let us introduce the notation 

OO OO OO OO C30 

Jm (a) = j dui J du2 ■ J dUm j dv2- ■ J dVm (143) 

— OO Ui Ujn-i 

m 

X e-^"? Y[ A (ui, ,,2 + U2, . . . ,i^m + u^) • 

1=2 

for the multiple integral entering into the expression of (a). Then we have 

frn («) = I — J™ (a) . (144) 

{2tt) 2 2! •• -(771-2)! 

The form of Jm (a) is reminiscent of the multiple integrals which appear in the 
theory of Gaussian random matrix ensembles. The following three lemmas show 
that in the three limiting cases Jm (a) can be explicitly evaluated in the form 
of Mehta integrals. 

Lemma 2 

^.(^-1) , , (-1)"'""" 2!---(m-2) , , 

hm a 2 Jm[a)^Im,i/2 j (145) 

where Im,i/2 the Mehta integral defined in i98\) . 

Proof Let us make a variable change under the integral ([25]) introducing new 
integration variables 

ipi = Ml, (146) 
ifi = Vi + Ui,i ^ 2, . . . ,m, (147) 
fii-i — avi,i = 2, . . . ,m. (148) 



In the new variables the integral ([25|1 can be written as 

^ m-l °f 

Jm {a) = m(rn-l)/2 11 / C^A^^Mr ^ 6^''' ^(^i , . . . , i ; «) , (149) 



where we introduce the notation 
g{^ii, . . .,Hm-i;a) 



dipi J dip2 J dip3 (150) 
J d(^™e-K^™+-+^?)A((pi,...,(p™). (151) 
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The function g{fii, . . . , ^m-i! ck) is bounded uniformly in a G R. 

\g{fll, ■ ■ .,flm-i;a)\ < I,n,l/2, (152) 

which can be shown by replacing the Vandermonde determinant under the in- 
tegral by its absolute value and extending the lower integration limits to minus 
infinity. Therefore the function under the integral in (|149p is uniformly bounded 
and integrable. By the dominating convergence theorem one can interchange 
the limit a — > oo and integration. Then, for the function 5(/i2, • ■ • ,A''m;Q!) we 
have 

hm g(/i2, . . . , Ai„; a) = 1^.1/2 j • (153) 

Remarkably the limiting value does not depend on the variables 
{/ii, . . . , /im-i}- Therefore the integration in (|149p can be performed indepen- 
dently for each i = 2, . . . ,m, each resulting in {i — 1)!. This yields p45p . 

Lemma 3 

Jm (0) = ; . Irn,l (154) 

(to — 1)!to! 
Proof Let us make the variable change 
Xi = ui 

Xt = -f Mj, z = 2, . . . ,TO. (155) 
Then the integral takes the form 



OO C30 



Jtn (0) = J J ^'^^ J du3 - J dUrn J (156) 

— OO Xl "^2 liTTi-l U2 

OO 

/m 
dXr«e-^Xi \[ {x^ - u^y-^ e-5X? A (xi, ■ • ■ , Xm) . 

The integrals over Ui, for i — 1, . . . ,m, can be evaluated step by step. First, for 
i = m, we have 



dUm I dXme '^^"^ iXm - Urn)"" ^ A (xi , • . • , Xm) (157) 



1 



TO — 1 



dXme ixm - u^-i)"' ^A(xi,...,Xm), 



which is done by changing the integration order. In the next step, the integral 
over Um-i can be calculated by parts: 



J dUm-l j dXm-1 J dXmS' 

1^771-2 Um-i Um-1 



2 A.m — 1 2 ■^m 
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X (Xm-1 - Wm-l)™ ^ (Xm - Mm-l)'" ^ A (xi, • • • , Xm) 

OO 

dXrn-1 [ dxr„e-^>^--i-^>^- (158) 



m- 2 

X (x A(x2,---,Xm) 



1) y" du^-i rfXm-l j dxme ^^^-1 



X (x m — 1 ^m— 1 m ^m— 1 )" A(xi,...,Xm)]- 

The second term cancels because of the antisymmetry of the Vandermonde de- 
terminant with respect to interchange of Xm and Xm-i- Iterating this procedure 
we remove (m — 1) integrals in the variables U2,. ■ . , Um ■ 

OO OO OO 

-OO XI XI 



xe ^'YliXi-XiT 2XiA(xi,...,Xm) 

i=2 



A symmetrization of the expression under the integral in the variables Xij i 
2, . . . , m, yields another Vandermonde determinant. Thus 



i(m — l) OO OO 



Jm (0) = [,^^ . j dxi j dX2 (160) 



-OO XI 



{{m-iy.y 

■■■J dx™e-Kx?+-+x^)|A(xi,...,x„.)|'. 



XI 



Finally we add this integral to the (m — 1) similar ones, obtained by interchang- 
ing xi with each of X2, • • • , Xm, and divide the sum by m. 



Jm{0) = \ , ' , / dXl (161) 

(m — l)!m! J 

— OO 

OO 

■■■ I dXme-Kx?+-+x^)|A(xi,...,Xm)r 

— OO 

This gives us the stated result. 
Lemma 4 

lim e-«=™(-i)/V„ (a) = J^-i,i ^"^^"" ^m-^ (162) 

((m-1)!) 
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Proof We start from the integral in ([^5]) and make a change of variables as 
follows, 

Xi = Mi + a + ui, i = l,...,m — 1 (163) 
Si = \a\{ui-ui), i = l,...,m-l (164) 
si = ui — a (to — 1) , (165) 

which yields the following integral expression, 

2tti(?ti — 1) ^ OO OO 

Jm (a) = / j ds2e-'^ j ds^e-'^ ■ ■ ■ (166) 



dsniC / dx2 ■ ■ ■ I dxm W {xi — Si — amY 

i=2 



.He 

i=2 



— GO 





OO 


OO 


J dx2 • • • 


/ 


-^am si- 


f am 












Sl 



52 Sm 
I I ■ ■ • I Xm + . . 

a \a\ 



Due to the presence of the Gaussian and exponential terms, the main contri- 
bution to the integral comes from finite values of si, . . . , Sm and a;2, . . . , Xm- 
Therefore, up to corrections of order of 0(1/ |a|), we can neglect the terms 
divided by |a|, and extend the lower limits of integration over X2, ■ ■ ■ ,Xm to 
— OO. The the integrals over S2, . . . , Sm decouple, and we can evaluate them to 
1/ (m — 1)!. The Vandermonde determinant becomes antisymmetric with re- 
spect to permutations of the variables X2 , . . . , Xm ■ As the integration is over the 
symmetric domain, we can leave only the antisymmetric part of the rest of the 
expression. The product nt=2 {^i ^ *i ^ then results in 

(-l)"^^^^-f^. (167) 
[m — 1)! 

After collecting the leading order terms from the first argument of 

A (si -I- am, X2 - ■■ , x,n) — {am)"'~^ A{x2 - ■■, Xm) , (168) 
the integral over si decouples as well, and yields y/n. We finally obtain 

. . N e 2 V27r(-1) 2 
Jm (ck) ~ — — 1 T (am) 

dX2--- / dXmd ^ ' \A{X2, ■ ■ ■ ,Xm)\^ ■ (169) 



Using the definition of the Mehta integral, ((98|) we obtain (|162p . 

The above lemmas, the formula for the Mehta integral and the definition 
(fTii]) of fra (a) estabhsh the results 
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6 Discussion of the results and conclusion 



The first result obtained in this paper is an expression for the survival prob- 
ability for m walkers in the SVW model. At the fixed parameter k < 1 and 
t ^ oo the leading asymptotics coincides with that for the usual 

VW [11]. This result is intuitively clear. In the case of the VW it is obtained by 
considering the evolution of independent noninteracting particles and reducing 
the number of its outcomes by dropping those realizations where the crossings 
of particle time space trajectories occur. Then, the asymptotics of the survival 
probability for two walkers follows directly from the diffusion law and the 

method of images. For m walkers, the method of images involves m(m — l)/2 
reflecting wall planes. Each wall brings with it a factor t~^/^ providing to- 
tal contribution f-™("i-i)/4 Pqj- ^ny fixed k < 1 the events of crossing occur 
with a finite, though changed, probability, so that this argument still holds. 
The survival probability changes by a constant factor dependent on k, leaving 
Fisher's law unchanged. As k approaches one this factor diverges, which in- 
dicates a qualitative change of behaviour of the survival probability. As the 
crossings become less probable, it finally saturates to the TASEP normaliza- 
tion constant. The transition between the two regions takes place on the scale 
(1 — k)^/! — const, where the effect of the diffusive spreading of particles be- 
comes comparable to the one caused by SVW interaction. In this case, the 
survival probability is expressed by the scaling function /,„(a), which has the 
SVW and TASEP asymptotics as limiting cases. 

The formulas for the survival probability in SVW can be reinterpreted in 
terms of the moment generating function of the time integrated particle current 
Yt in TASEP, which, in turn, can be used to construct the distribution of the 
same quantity. The cases of generic k < 1 and k > 1 correspond to the tails 
of the probability distribution of Yt at the large deviation scale, characterizing 
positive, (Yt — {Yt))/t > 0, and negative, {Yt — {Yt))/t < 0, deviations respec- 
tively. The positive tail has the form specific for the current distribution for m 
free independent Bernoulli particles, while the form of the negative tail looks 
like that of the distribution for one Bernoulli particle that makes m steps at 
a time or m particles jumping one step synchronously. Such asymmetry re- 
veals different mechanisms of positive and negative fluctuations. For positive 
deviations m particles have to be " accelerated" independently. The main contri- 
bution to the negative deviations comes from the events when the first particle 
decelerates all particles following behind. Our results are to be compared to 
the ones obtained by Derrida and Lebowitz [S] (see also [7]) for the large devi- 
ation function of the TASEP current on a ring. In particular, they have shown 
that for large m the non-universal tails of the current probability have the form 
P{Yt/t = w) ~ e"^+ and P{Yt/t = w) - e^-^^/™)' for positive and neg- 
ative deviations respectively. Our results have the same functional form even 
for finite m with H^{v) = H+iv) being a simple large deviation function of 
the Bernoulli process. Remarkably, such a mechanism survives on the infinite 
lattice, despite the particles drifting apart from each other in the course of time, 
so that they meet less and less often. The acceleration-deceleration asymmetry 
was also observed in the large deviations of the distance travelled by individ- 
ual particles in TASEP studied in for a general case of particle dependent 
hopping rates. There the negative large deviations do not depend on the order 
number of a particle whereas the positive ones do. 



33 



The result obtained for the SVW in the transition region (1 — k)^/! = const 
provides us with the hmiting distribution of the particle current measured at 
the diffusive scale, \Yt — (Yt)! ^ ^/t. The distribution obtained is parameter 
free, dependent only on the number of particles. We expect that it is a uni- 
versal distribution for the systems of particles performing a driven diffusion, 
independent of the details of microscopic dynamics. The distribution has a 
skew, non-Gaussian form, with tails matching the large deviation behaviour 
asymptotically. 

Several directions for future work can be mentioned. First direct contin- 
uation of the present paper is an asymptotic study of the function {a) as 
m ^ oo. Considering a certain common scaling with the variable a is expected 
to give a new universal scaling function characterizing the KPZ class. Note 
that in the limiting cases this function reduces to Mehta integrals, which are 
the probability normalization constants for Gaussian random matrix ensembles, 
namely the Unitary and Orthogonal ensembles. The asymptotic evaluation of 
these integrals performed in the random matrix theory resulted in non-trivial 
densities of critical points, from which comes the leading contribution to the 
integrals [32] . It would be interesting to see how these densities transform from 
one to another as the parameters vary between the limiting cases corresponding 
to different ensembles. It would also be interesting to study the TASEP con- 
fined in a ring. The starting point of this analysis could be the recently obtained 
expressions for the Green functions [35]. In this way one could obtain a scaling 
function that characterizes the behaviour of KPZ interfaces in finite systems. 

Another possible development of the present article is a generalization of the 
above mentioned results for the probability distributions of the distance trav- 
elled by an individual particle in the TASEP and the corresponding correlation 
functions to SVW case. Note that similar results exist also for the VW 1 . In 
both cases the appropriately rescaled distribution of the distance travelled by 
a single particle starting from a half filled lattice is shown to converge, to the 
so-called Tracy- Widom distributions j43j, which appear in the random matrix 
theory as a distribution of maximal eigenvalue in the Gaussian ensembles, uni- 
tary in case of TASEP and orthogonal for VW. The SVW model establishes 
a bridge between these two cases. However, its Green function has neither a 
Toplitz form like in VW nor a special structure like in TASEP, which allowed 
Sasamoto, [3S], to reinterpret it again as a problem of the VW and finally to 
represent its evolution as a determinantal point process. Therefore, a significant 
extension of the existing techniques is in order. In this connection we should 
mention the recent advance for the Partially Asymmetric Simple Exclusion Pro- 
cess ^43]- [35], which was made only on the basis of the Bethe Ansatz solution 
without use any free fermionic representation like VW. 

An interesting example of SVW has been proposed recently by Johansson 
[22j in his analysis of a domino tilling problem on the Aztec diamond known 
as the arctic circle problem. It was shown that the domino configurations are 
in one-to-one correspondence with trajectories of an n-particle process which 
is defined as follows. At each discrete time step a particle jumps forward with 
probability q or stays put with probability p = 1 — q. If the next site is occupied, 
the probability to stay put is 1 — q{l — k) as in the SVW model. In addition, 
after each step, a particle i can be translated back to the distance Si with 
probability q^^ provided that si < Xi — Xi^i for all i. If one chooses k — ~q, 
the model belongs to the free fermion class and its transition probabilities admit 
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a determinant representation. It has been shown in [22 that the position of the 
first particle is described by the Airy process in the thermodynamic hmit for the 
domain wall boundary conditions in the domino tiling problem. By similarity of 
the models, one may expect that the extremal statistics of the SVW model also 
exhibits a kind of Tracy- Widom distribution for appropriate initial conditions. 
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